Astronomy & Astrophysics manuscript no. great'atm 
February 25, 2013 



© ESO 2013 



GREAT/SOFIA atmospheric calibration 

Xin Guan^, Jiirgen Stutzki^, Urs U. Graf^, Rolf Glisten'^, Yoko Okada^, Miguel Angel Requena Torres^, 

Robert Simon^, Helmut Wiesemeyer^, 



(N 

o 



l> 



^ I. Physikalisches Institut, Universitat zu Koln, Ziilpicher Str. 77, Koln, Germany, 50937 
^ Max-Plank-Institut fiir Radioastronomie, Auf dem Hiigel 69, Bonn, Germany, 53121 

To appear in the A&A special issue for GREAT/SOFIA 

ABSTRACT 

Tiie GREAT observations need frequency-selective calibration across the passband for the residual atmospheric opacity 
at flight altitude. At these altitudes the atmospheric opacity has both narrow and broad spectral features. To determine 
the atmospheric transmission at high spectral resolution, GREAT compares the observed atmospheric emission with 
atmospheric model predictions, and therefore depends on the validity of the atmospheric models. We discusse the 
problems identified in this comparison with respect to the observed data and the models, and describe the strategy 
used to calibrate the science data from GREAT/SOFIA during the first observing periods. 
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1. Introduction 

The German REceiver for Astronomy a t Ter ahertz fre- 
quencies (GREATfl (iHevminck et all . I2012D onboard 
the Strat ospheric Observa t ory for Infrared Astron omy 
(SOFIA) (|Becklin fc Gehrzl . [2009t lYoung et all . [20121) ob- 
serves at typical flight altitudes of 8 km (where telescope- 
and mstrument set-up starts) up to 14 km. The atmospheric 
transparency at these altitudes is high (though not unity) 
except near strong absorption lines of species such as wa- 
ter and other trace constituents (O3, CO, N2O, etc.). The 
typical amount of water vapor in the atmosphere above the 
observatory (precipitable water vapor, pwv) is in the range 
of a few /im at very good observing conditions and up to 
50 ^m and higher at bad conditions. The atmospheric pres- 
sure, rapidly decreasing with increasing altitude, results 
in narrower absorption lines with increasing altitude. As 
a result, the atmospheric opacity at flight altitude shows 
many narrow spectral features, typically with a width of 
about 100 MHz down to a few MHz (varying with species 
and altitude at which the particular species absorbs). In 
addition, extended line wings from very strong absorption 
lines, in particular of water, and collision-induced absorp- 
tion (CIA), where the collisions of the absorbing species 
with other molecules generate a temporary dipole moment, 
contribute to a spectrally smooth broad-band opacity. 

The correction of airborne astronomical data for the at- 
mospheric transmission is thus much more challenging than 
for ground-based observations, because of the small-scale 
variations resulting from the narrow atmospheric lines. 
Direct measurement of the frequency-dependent transmis- 
sion is in principle possible, but is limited by the inherently 
higher noise in narrow frequency-bins. Hence, the calibra- 
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tion has to rely on appropriate models of the frequency- 
dependent atmospheric opacity and the proper approach is 
to determine the transmission from model fits to the ob- 
served sky emission. This approach has been implemented 
and successfully used in recent years a t ground-based high- 
altitude observato ries such as APEX (iGiisten et al.l . l2006[) 
and NANTEN2 (jKawamura et al.l . |2005D in Chile. The 
proper description of the frequency-dependent atmospheric 
transmission is particularly important for absorption mea- 
surements against continuum background sources, because 
the background emission is also modulated by the atmo- 
spheric transmission, which can therefore interfere with the 
absorption profiles from the astronomical sources. 

We first discuss the calibration approach in section [5] 
and its implementation for GREAT on SOFIA in section[3l 
GREAT allows us for the first time to simultaneously ob- 
serve two widely separated frequency bands in the THz- 
regime at high spectral resolution. Accordingly, it has come 
as no surprise that the atmospheric models currently in use 
for predicting the THz-regime emission and absorption of 
the atmosphere turn out to have clear deficiencies in prop- 
erly modeling the atmosphere above SOFIA to the degree 
necessary to consistently calibrate the GREAT data. These 
problems are discussed in section HI We conclude by de- 
scribing the pragmatic approach chosen for the calibration 
of the GREAT data of the first observing season and the 
perspectives for future improvements in section [5j 



2. Approach to atmospheric calibration 

Following the notation of Appendix |A1 the sky brightness 
on an equivalent antenna temperature scale along a line-of- 
sight at elevation El for a plane parallel atmosphere, i.e. 
taking for the airmass the approximation A = l/sin(£'/), 
is given by 



T, 



AjSky,^ 



Ju {T{tI)) X 
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(1) 
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with the atmospheric transmission 

ta,u = exp {-Tz^vl SVCl{EI)) , 



(2) 



(3) 



the frequency-dependent zenith opacity tz,v^ and the defi- 
nition of the effective sky brightness 



^shy.y — ^A^shy.vj \^ ^aM} 



(4) 



according to equ. IA.7I 

An often used approximation is that the emission occurs 
at the ambient temperature of the atmosphere near the 
observatory, Tamt, so that 



-^sky.u ^^ ui/\-^amb)j 



(5) 



ahhough this assumption may not be particularly well justi- 
fied for absorption lines from the higher atmosphere, which 
is warmer than at flight altitude (see Section [3]). 

In principle, the above equation can be inverted to cal- 
culate the sky transmission at each frequency from the mea- 
sured sky brightness. However, because the atmospheric 
calibration has to be frequency specific, deriving the opac- 
ity from the sky-emission in each frequency bin results in 
quite large errors due to the limited signal-to-noise ratio 
in their narrow band-width. In addition, a heterodyne in- 
strument receiving double-sideband (such as GREAT) folds 
the signal- and image-sideband on top of each other, so 
that the emission from the two sidebands cannot be sepa- 
rated. To overcome these problems, one has to (i) use an 
atmospheric model to predict the atmospheric opacity as a 
function of frequency at the given observing altitude and 
line-of-sight, parameterized by a few dominant atmospheric 
parameters, such as the precipitable water vapor and the 
ambient temperature, (ii) derive these parameters from the 
comparison between the observed emission and the model, 
and (iii) then use the model opacities, calculated with these 
best-matching parameters, to correct for the attenuation of 
the astronomical signal. 

This can be done either for individual spectrome- 
ter channels to only de-convolve the image- and signal- 
sideband opacity, or, to increase the signal-to-noise ratio, 
over selected sub-bands or the complete reception band- 
width of the receiver (which in the case of GREAT might 
include the two simultaneously observed THz-bands LI and 
L2) . A necessary condition for the consistency of the cali- 
bration approach is that the atmospheric parameter(s) de- 
rived for individual channels or sub-bands are identical to 
within the noise and/or that the residual of the compari- 
son between the atmospheric model and the observed sky- 
brightness is within the observed noise level. If this is not 
met, the atmospheric model obviously does not describe the 
actual atmosphere satisfactorily (or the measurements are 
corrupt). 

The zenith opacity above the observer at altitude h de- 
pends on the altitude profile of the abundances of the at- 
mospheric species that are responsible for the absorption, 
which is driven by the complex chemical processes in the at- 
mosphere, as well as by the pressure profile (influencing the 
pressure broadening) and temperature profile (determining 
the population of the relevant species within the absorbing 
states). All these may vary as a function of geographical 
position and local weather parameters. To arrive at a man- 
ageable problem, one has to make additional assumptions. 



such as that the altitude variations of abundances, temper- 
atures, and p ressure arc given by a standard atmosphere, 
e.g., the lU.S . Standard Atmosphere (1976), which specifies 
the altitude variation on a pressure altitude vertical scale 
(in the following, altitude often refers to pressure altitude). 
Similarly, it specifies the total abundances and abundance 
variation with altitude of the atmospheric species of inter- 
est, except for a few, where their distribution has to be 
adapted to the specific observing conditions, as discussed 
below. In particular the column of water molecules above 
the observer, pwv, may change rapidly, because it is largely 
determined by local weather. Hence, a useful parameteriza- 
tion of the atmospheric absorption is given by splitting the 
wet, i.e., water absorption, and the dry part as 



To,,. = b„{Pambl Tamb, X-j,.. .]) X pwV + 



(6) 



where the wet and dry coefficients are a function of the 
ambient pressure Pamb (or altitude), and additional atmo- 
spheric parameters such as the ambient temperature Tamb, 
abundances Xj, and possibly other parameters that, e.g., 
specify the altitude profile of these abundances. A par- 
ticular atmospheric model can then be represented by a 
multi-dimensional look-up-table of the h, c -coefficients vs. 
frequency, (pressure) altitude, and additional parameters, 
which are derived from the various atmospheric models by 
fitting the model opacities, calculated at two values oi pwv 
and at different frequencies and altitudes. The details on 
deriving the b- and c-coefficients are given in Appendix |^ 

2.1. Atmospheric models 

Several atmospheric models are in use within the radio- and 
FIR-astronomy community, both for ground-based observa- 
tories and air-borne applications. For the GREAT/ SOFIA 
calib ration we have i nvest igated the use of AM (iPaine . 
"mit), ATRAN (|Lordl . fT99^. and MOLIERE ([Urban et al. . 
2004). The ATM-model (jPardo et aU l200l|), in use for 
ALMA and other ground-based observatories, in its pub- 
lically available form does not include the higher excitation 
water lines and hence is of limited use for frequencies above 
1.5 THz. 

Figure (TJ comparing the wet- and dry- atmospheric 
coefficients of the three models, shows that ATRAN at 
present does not include the collision-induced absorption 
(CIA) from N2 and O2, which is responsible for the quasi- 
continuum opacity shown in the other two models. The at- 
mospheric opacity due to these CIA-processes differs by a 
factor of w 2 between MOLIERE and AM, being higher 
for MOLIERE. All three models agree on the wet atmo- 
spheric opacity, although ATRAN has some small errors 
in the frequencies of a few higher excitation water lines. 
Because ATRAN does not include the dry-atmosphere con- 
tinuum opacity, and the dry-atmosphere contribution to the 
opacity due to CIA-processes is already too high compared 
to the observed sky emission (see below) even for the AM 
model, we used the AM model for the standard GREAT 
calibration. 

To calculate the atmospheric transmission for differ- 
ent pwv-vahies, we modeled the water mixing ratio as 
an exponential function with altitude in the high tropo- 
sphere and low stratosphere, a nd with a constant , but 
low value higher above, following iMastenbroold (jl968t) and 
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Fig. 1. Comparison of the dry (top) and wet (bottom) opacity coefficients as defined in eau. 151 [A. 291 and IA.27I for the 
frequency range of GREAT channels LI (left) and L2 (right), calculated from the atmospheric models AM, ATRAN 
and MOLIERE (see text), showing the opacity for a pressure altitude of Pamb — 133Torr and a pwv — 10/im. The 
largest discrepancy is in the quasi-continuous dry-atmospheric opacity due to collision- induced absorption (CIA), which 
is completely neglected in ATRAN and differs by about a factor of 2 between MOLIERE and AM. 



iHerbin et al.l (|2009f ) . The detailed treatment is presented in 
Appendix \^ 

The nitrous oxide, N2O, pr ofile is known from 
ground- and balloo n observa tions ([Emmons et al.l . 119941 : 



20081^ 



We adapted 

uses a N2O 

of 60 hPa, 



iCamv-Pevret et all . 12004: Strong et al 
a simplified profile for use with AM. It 
mixing ratio of 0.3 ppm below an altitude 
and 0.1 ppm above. This is a rough approximation of 
many measured profiles, s uch as that given in Figure 1 of 
ICamv-Pevret et all (|2004D . N2O does not produce particu- 
larly strong atmospheric features, but we adapted an N2O 
profile based on real observations, to compare AM with 
ATRAN, which produces far too strong N2O features com- 
pared with the observed sky emission. 

The ozone, O3, profile changes rapidly with geographic 
positions and with season. The profile actually used in 
the GREAT calibration is obtained by fitting sky mea- 
surements with the NANTEN2 telescope in Chile to 
properly reproduce the strength of the ozone lines. The 
GREAT/SOFIA calibration shows that the measured ozone 
lines are adequately reproduced by the model, so that 
this altitude profile is also applicable here (although a few 
SOFIA measurements show an increased ozone abundance 
compared to this model, which does not affect the overall 
calibration, because the ozone lines are very narrow spectral 
features, but has to be taken into account for absorption 
measurements against a continuum source). 



3. Implementation for GREAT 

The calibration of atmospheric transmission for 
GREAT/SOFIA is implemented in the task kalihrate 
as part of the kosmasoftware observing software package. 
The same scheme is used for the submm-wave data at the 
NANTEN2 observatory in Chile. It fits the observed, cali- 
brated differential sky-hot spectrum with the atmospheric 
model brightness as given by eaus. lA.48] and [S.47l fand uses 
the approximation uj = fl = and Tsky,w — Ju(Tsky) as 
detailed in Appendix [2| ■ The best- fit value is thus deter- 



mined with high precision, making use of the full spectral 
information observed, i.e., with a good signal-to-noise 
ratio. The user can control whether the fit is performed 
as a "common" fit to all simultaneously observed spectra 
(several receiver pixels at several frequency bands) or 
individually for each spectrometer. The atmospheric 
transmission for this best-fit-value is then applied to the 
full spectra channel by channel, thus keeping the full 
frequency resolution of the atmospheric features (to the 
resolution with which the frequencies where sampled when 
generating the 6-,c-coefficient look-up tables). 



The approximation Tsky,v — J^u{Tsky) assumes that 
all opacity originates at the same temperature. To check 
this assumption. Figure [5] shows the effective sky brightness 
Tsky,v (equ. 131), calculated with the AM model from the 
model brightness and opacity for a frequency range near 
a strong water line at two different pww-values. The dry- 
atmospheric lines (shown in the middle panel in Figure H]) , 
as well as the narrow core of the high excitation water line 
(see lower panel in Figure [2]), originate in the higher strato- 
sphere, where the temperature increases well above the am- 
bient temperature at the observer's altitude in the lower 
stratosphere. At the low pwv of 3 /im, the effective sky 
brightness temperature in these lines is therefore substan- 
tially higher than the Rayleigh-Jeans-(R-J) corrected am- 
bient sky temperature of about 172 K. With higher pwv, 
the increasing opacity at altitudes, and hence tempera- 
tures, near the observer's altitude absorbs these emission 
lines from the high atmosphere, so that the excess effective 
sky temperature near these lines is reduced and the effec- 
tive sky temperature approaches the ambient temperature. 
The excess is typically up to about 10 K at very low pwv. 
Neglecting it results in a slight underestimation of the sky 
brightness in the cores of these lines for the model fit to 
the observed brightness. Because this affects only narrow 
frequency ranges, it has a negligible effect on the derived 
fit parameters. 
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Fig. 2. Effective sky brightness, Tgky, (top panel) in a rep- 
resentative frequency window near a water line and several 
dry-atmospheric lines originating in the higher, warmer at- 
mosphere, calculated with the AM model at a pressure al- 
titude of 133 Torr and an ambient temperature of 200 K, 
for two values of pwv, 3 (red) and 30 /xm (blue). The 
lower two panels show the dry atmospheric opacity (i.e., 
the c-coefficient; middle) and the wet atmospheric opacity 
at 10 /im pwv, i.e., b x 10 /im; bottom), similar to Fig. [T] 
to identify the spectral features causing the excess effective 
brightness temperature. 
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Fig. 3. Example Tsky - T^hot spectra from simultaneous 
observations in the LI (red) and L2 (green) band, showing 
two ozone lines (the broad features toward the middle of 
the LI and L2 passband) and two CO lines (the narrow 
peaks on the left). Note that the two CO lines appear at 
the same velocity on a vlsr scale, because the receiver was 
tuned to two CO lines for the astronomical observations. 
The two black curves are independent fits of Tsfc^ - That 
for the two bands. The red and green lines between the 
measured spectra are the result of a common fit, which 
simultaneously fits a single value of pwv to both bands. In 
a common fit, a consistent value oi pwv cannot be derived 
for both bands. 
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Fig. 4. Time series of hot-cold load-calibration scans (in 
the GREAT L2 channel), displayed across the central part 
of the IF band of width 1.4 GHz, calibrated against the pre- 
vious load scan and normalized to the R-J-corrected hot- 
cold brightness difference. It shows that the receiver gain 
drifts by up to a few percent over time scales of several min- 
utes. The sequence number and time offset to the previous 
scan is given in the bottom left corner of each panel. 



4. Application to GREAT data 

This calibration scheme, when applied to the GREAT data, 
is not able to self-consistently fit the atmospheric emission 
simultaneously for the LI and L2 bands, whereas indepen- 
dent fits to the individual receiver bands converge well, but 
on solutions with different values of pwv for each band. 
Comparison of the fits with the observed sky-hot spec- 
tra shows that the narrow spectral features resulting from 
the dry atmosphere and the line features of water are re- 
produced well. This indicates that the problem lies in the 
quasi-continuum opacity of the dry atmosphere, i.e., the 
approximately constant offset of the c-cocfficients visible in 
Fig. [T] (resulting largely from the CIA-processes of N2 and 
O2) predicted in the models. 

A closer look confirms this and shows that this quasi- 
continuum opacity of the dry atmosphere is too high. 
Indeed, independent fits to the emission in the two re- 
ceiver bands for observations with very good transmission 
converge to a formal fit result with slightly negative pwv- 
values: because the model opacity cannot be reduced below 
the dry-atmospheric contribution, i.e., the (too high) value 
of the c-coefficient, this excess opacity has to be compen- 
sated by negative pwv, and hence a negative wet opacity, 
b X pwv, to match the observed low brightness of the sky 
at these low opacities (because negative values for pwv are 
unphysical, the fit result is reset to pwv = in these cases) . 

Observations at higher overall opacity, i.e., low elevation 
and higher pwv, show a similar discrepancy. Figure [3] shows 
such an example for a tuning that observes the CO (11- 
10) line in LI and the CO (16-15) line in L2. Separate 
fits (black lines in Figure [3]) give pwv=12.3 /im for the LI 
band and 35 /im for the L2 band. The excess opacity in the 
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model needs to be compensated in the fitting by a more 
reduced wet opacity, hence reduced pwv values, than would 
be needed with the proper dry opacity. The scaling factor 
between the wet opacity and pwv, namely the 5-coefficient 
as defined above, is typically higher by a factor of 2 to 3 for 
the L2 band, so that the excess dry opacities in the model 
result in accordingly scaled, i.e., smaller differences in pwv 
for the fitted pwv-vahie for the L2 band compared to LI. 

Similarly, a common fit to both bands, now forced to 
using a unique value for pwv, results in a pwv-vahie some- 
where in between the two discrepant values fitted in the 
separate fits. It overestimates the opacity in the band with 
lower opacity and underestimates the opacity in the other, 
higher opacity, band as shown by the green and red fit line 
in Figure[31 which corresponds to a best-fit value of 21.4 fim 
pwv. 

The above holds for tunings where the wet atmosphere 
does not contribute any line feature, but merely a contin- 
uum opacity (with a slight slope). In the rare cases where 
the wet atmosphere also contributes with a line feature, 
even the single frequency band fits may not converge to a 
consistent result: the line intensity forces a certain amount 
of water and hence also continuous opacity, so that the to- 
tal opacity, including the dry opacity overestimated by the 
atmospheric model, overshoots the range that is consistent 
with the observed atmospheric emission. 

To confirm these findings, we modified the fitting pro- 
cedure to include additional fit parameters oli and ol2, 
namely a frequency dependent dry-offset, so that equ. [S] 
changes to 



Tu = bt, pwv + Cy - 



OLI: 
0L2: 



V in LI 

V in L2 



(7) 



This gives a more consistent fit, but not with unique offsets 
Oy for different tunings and an ambiguous separation into 
pwv offsets and values for o^. Similarly, one could use an 
empirical atmospheric model such as a linear combination 
of two of the models above. In practice the two extra pa- 
rameters of course give a better fit to the data, but allow 
different combinations of their values and thus do not con- 
strain the fit toward a better insight into the atmospheric 
model deficiency. 

In principle, one can estimate the value by which the 
model-predicted quasi-continuum dry opacity is too high 
compared to the observations by analyzing the data taken 
at very low pwv, i.e., where the formal fit gives negative 
vahies for pwv. However, this requires a very precise mea- 
surement of the low residual sky-brightness at these low 
opacities. Despite the excellent stability of GREAT, doc- 
umented by the Allan- va r iance minimum times of around 
100 sec (jHevminck et al.l . l2012l ). slight drifts in the total 
power gain profile caused by, e.g., mechanical deformation 
under gravity of the diplexer optics or residual drifts with 
temperature in the receiver electronics, limit the precision 
with which these measurements can be made. Figure |4] 
shows that these gain variations on time-scales of several 
minutes, i.e., between subsequent hot-cold gain-calibration 
scans, are on the order of up to a few percent. These rel- 
ative errors in the gain translate into absolute errors in 
the derived opacities, because the corresponding sky-hot- 
observations are affected at the same level. Note that ab- 
solute errors in the opacity of this order are perfectly tol- 
erable for the calibration of the atmospheric transmission. 



because they contribute a small relative error of similar 
magnitude to the transmission, and hence to the inten- 
sity of the calibrated spectrum. But they compromise a 
more precise analysis of the atmospheric model discrepan- 
cies discussed above, because we are talking about small 
offsets, i.e., corrections to the c-coefficient on the order of 
a few times 0.01. Only the few sky- measurements immedi- 
ately following a load-measurement avoid this problem, but 
there are too few at very low opacity within the data set 
from the SOFIA basic-science flights, and these are scat- 
tered over different receiver tunings, so that they do not 
give sufficient statistics to derive a consistent estimate of 
the model discrepancy. 



5. Summary and outlook 

Given the problems in the atmospheric models as discussed 
above, the GREAT data from the basic science series of 
flights were calibrated using individual fits to each receiver 
band. Although these give inconsistent values for the best- 
fit pwv for simultaneous observations, the relevant entity 
entering the calibration of the data is the transmission and 
not pwv. The calibration thus does not take advantage of 
the higher signal-to-noise that would be achievable with a 
common fit to the full reception bandwidth of the receiver, 
but it is nevertheless sufficiently precise to not introduce 
any significant errors (depending on the detailed circum- 
stances of the observing conditions, we estimate the error in 
the atmospheric transmission thus derived to be on the or- 
der of a few up to a maximum of about 10%). Correcting for 
the atmospheric transmission by applying the frequency- 
dependent opacity as given by the atmospheric model, in- 
cluding the narrow spectral features resulting from the dry 
atmosphere at full spectral resolution, is crucial, however, 
to remove residuals from these features in the resulting as- 
tronomical spectrum, in particular for absorption measure- 
ments against a continuum source. 

To improve the calibration, the inconsistencies in the 
atmospheric models with regard to the level of broad- 
band, quasi-continuum opacity that we identify with the 
GREAT/SOFIA observations in the THz-regime have to 
be overcome by developing a better understanding of 
the broad-band, quasi-continuous opacity contributed by 
the CIA-processes; in particular, its broad-band variation 
across the frequency range covered by the different channels 
of GREAT needs to be properly understood. In parallel, a 
careful analysis of the sky- measurements available from the 
basic-science flight series, but also dedicated future obser- 
vations tuned to particular frequency settings with strong, 
but narrow water lines from the slightly higher atmosphere 
that may allow us to independently fix the pwv value and 
thus to resolve the ambiguity between dry continuous and 
wet continuous opacity contributions, will help to quantify 
and possibly correct the discrepancy in the models. 
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Appendix A: Atmospheric calibration of heterodyne astronomical observations 

(to appear online only) 

This appendix is a complete compilation of the calibration scheme for astronomical observations with a heterodyne 
receiver, taking into account the special problems that appear at highest frequencies, i.e., at the edge or beyond the 
validity of the common Ray leigh- Jeans (R-J) approximation, where the R-J correction across the receiver input band 
has to be taken into account with the high IF-bandwidths in use at present instruments. It also describes a general 
scheme for incorporating the correction for atmospheric transmission, which for observations at high altitudes from the 
ground, but in particular from air-borne or balloon-borne altitudes requires a detailed treatment of the spectrally narrow 
atmospheric features. 

The terminology follows, to the extent possible, the notation defined in lKutner fc Ulichl (|1981[ ) and iDownesI ()1989t ). 

A.l. Specific intensity and equivalent brigiitness temperature 

The specific intensity I^ is expressed as an antenna brightness temperature Ta,v m the usual radio-astronomical notation, 



For a black-body radiation field with a physical temperature T, 



B.{T) 



2hv^ 



^hv/kT 



the brightness temperature is thus 



■^'<^)-^«'<^)-T 



^hu/kT __ -^ 



(A.l) 



(A.2) 



(A.3) 



Note that expansion in |f < 1 gives J^{T) ~T-^[l- \^ -h O {{^f) ] . The difference between two black-bodies 
at temperatures T2 and Ti thus is 



MT2)-MTi)^{T2-Ti) 



I f hv\ f hv 



= {T2-T,)[l-6]. 



(A.4) 



With -^ = 48 K X t^ [THz] , the correction term for a frequency oi v = 2THz and temperatures of 300 K and 77 K hence 
is (5 = —0.066, small but still significant, whereas at 345 GHz it would only he 5 = —0.0019. 

Another quantity of interest is the ratio of black-body intensities between the image- and signal-sideband. With 
Vi — Vs :L I'lF, we obtain 



JuAT) '^ ^\ ZkT) kT '■ 



(A.5) 



which, for v — 2THz, vjp — 4GHz and a representative low temperature of 77 K gives an upper-/lower-sideband ratio of 
1.0035, i.e. within 0.35% of unity. 



A.2. Antenna temperatures 

We consider the general case of a heterodyne mixer that is sensitive in both sidebands and looks out to the sky, which 
emits an intensity equivalent to an antenna temperature T^^sfcy,!/, partially picks up intensity from ambient material with 
a physical temperature Tamb, and possibly has a sideband filter that is terminated at a physical temperature of Tterm- 
The mixer then sees in a narrow frequency band, centered at a signal sideband frequency Vg and an image sideband 
frequency Vi, a total intensity of 



Ta — Vmb ta.i^s ^sXs TmB.u^ + Vmb ia.Ui CriXi Tj^iB^u- + 

V-L Jamb) v-L ^a.Us /^S'^s -^sky^Us i V ^a,i>i j^i'^i -^sky^Ui ~r 

Jamb \y^ S'^ s <J u s \-^ arab } ^^ ^i^i <J Ui\-^ amb)\ '^ 

(1 ~ Xs)Gs Jvs (Tterm) + (1 ~ Xi)Gi J^^ (Tterm) ^ 

where we have introduced the following definitions: 

Jamb fraction of ambient material in sky beam 

•qmb niain beam efficiency 

Tmb,s,i source main beam brightness temperature in signal(image) sideband 

Tsky,u effective sky brightness temperature at frequency i/. 



(A.6) 



(A.7) 
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^sky^u • ^ A^sky-u j \^ ^ci.u) 

ia,i/ atmospheric transmission at frequency ly 
Gs,i signal(image) sideband gain: Gs + G^ = 1 
Xs^i signal(image) response in sky beam, 

with coupHng to sideband filter termination: (1 — Xs^i). 

Here, we have assumed that the coupling factor to ambient and the main beam efhciency are independent of the reception 
sideband. 

A. 3. Types of receivers 

We can distinguish the following two cases straightforwardly: 

— case A: a receiver with no sideband filter has 

Xs — Xi — 1 

T™ /^ 1/^ 1^/^ 1 r^ 

l+Gi/G 



=> (jsXs + C^iXi — (js + Gj — 1 => Gs — 1 — Gi — 1+G- / ~ 



GREAT on SOFIA falls into this category. 
— case B: a receiver with a sideband filter, i.e., non-zero coupling to the sideband termination, has 

Xs +Xi = I ^ Xs ^ I - X, = —— — - — . (A. 9) 

-L I Xi f X g 

Other hardware implementations correspond to different values of Gs,i and Xs,i- 

A. 4. Main-beam brightness temperature and forward beam brightness temperature 

In the above, we have expressed the source signal as a main-beam brigfitness tem,perature, i.e., as the convolution of the 
source brightness distribution on the sky with the main-beam profile of the telescope. The coupling of the detector to 
the telescope main beam is the main-beam efficiency rimt- This is appropriate if the telescope has a clean main beam, 
thus avoiding additional pick-up of source intensity in side-lobes (error-beam pickup) . This is also a good description if 
the source is compact so that it has only intensity within the main beam. 

Alternatively, the source signal can be expressed as a forward beam brightness temperature T^ j^, i.e., as the con- 
volution of the source brightness with the full forward antenna response. The relevant coupling efficiency is then the 
forward efficiency ry/md, which in the above nomenclature is rjf^d = 1 — famb- The term rjmb Tmb,v then is replaced by 

Vfwd T^^, = (1 — Jamb) T^^i,- 

A. 5. Sky brightness and opacity 

A. 5.1. Opacity 

The sky transmission along a line-of-sight at elevation El for a plane-parallel atmosphere, i.e., approximating the airmass 
hy A = 1/ sm{El) (the latter assumption being valid except at extremely low elevations), is 

fa,. -e-"-/"'"(^') (A.IO) 

is given by the zenith opacity Tj/, where we drop the index Z in the following. 

The zenith-opacity of a particular species ^ is the integral along the line-of-sight from the observers altitude h to 
space over the abundance of that species times its absorption coefficient. The abundance at altitude s is expressed as the 
volume-mixing ratio 

x^{s) ^n^{s)/n{s) (A. 11) 

relative to the total density n{s) (we closely follow the nomenclature of the AM handbook ()Painell2011f )V The general 
expression for the opacity of a particular atmospheric species ^ with number density n^{s) = a;j(s)n(s) is then the 
line-of-sight integral over the product of its abundance times its pressure-, P{s) , and temperature-, T(s), -dependent 
molecular absorption coefficient, 

/•OO 

T^,^(h)= x^{s)n{s)k^{i',P{s),T{s))ds. (A.12) 

Jh 

The total opacity is obtained by summation over all species, 

T,{h)^Y.^,^dh)- (A.13) 
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The details of the microphysics, i.e., the kind of absorption (single line transitions, collision-induced absorption, etc.) are 
contained in the molecular absorption coefficients. For collision-induced absorption (CIA) the coefficient is proportional 
to the the abundance of the collision partner(s) £,, i.e., 

k^ - x^{s)nis)^'^Yiiy,Pis),Tis)), (A.I4) 

which for self-broadening (^ — ^) leads to a quadratic dependence on the abundance of the species. 

Instead of altitude, the equivalent pressure can be used to specify the height in the atmosphere (pressure altitude) 
through the hydrostatic equation 

dP 

n(h)dh= , (A.15) 

rriair g 

with the mean molecular weight ruair of the atmosphere and the gravitational acceleration g. 

We will show below that the zenith opacity of the atmosphere above the observer can be separated into a wet and a 
dry part (where the dry part may contain opacity from the high water column in the stratosphere with constant volume 
mixing ratio that does not change with local humidity and daily weather). Thus, the opacity can be written as 

T^{h) = K{h)pwv + c^{h). (A.16) 

A. 5. 2. Sky brightness 

The sky brightness is the integral along the line-of-sight at elevation El over the source function times the attenuation, 
which for a plane-parallel layering, approximating the airmass as A = 1/ sm{El), reads 

TA^skv,u = / ' J{T{tI)) exp {~tI / sin El) dT'J sin El. (A.17) 



From this, the effective sky brightness as defined in eau. lA.Tl can be calculated with a given line-of-sight transmission taM 



as 



Tsky,v = TA,sky,vl (1 — ta,v) ■ (A. 18) 

A. 6. Parameterization of tlie atmospheric water vapor 

For water, a reasonable altitude distribution in the stratosphere and upper troposphere is given by a low, constant mixing 
ratio, domina nt in the upper stratosphere, plus an exp onential term, extending through the upper troposphere and lower 
stratosphere (jMastenbrookL Il968t iHerbin et al.l . 120091 ). We thus write 

xh2o{s) = Xc + Xe eyij>[-{s - he)/hs] (A.19) 

= Xc + x'^exp{-s/hs), 

with Xg = Xe exp(/ie//is). The typical value for Xc is 2.5 10~^ ([Mastenbrookl 119681 : iHerbin et al.l . I2OO90 : he or x'^ can be 
adjusted to match the desired column of water vapor above the observer. Eau. lA. 191 gives a total column density of water 
above altitude h of 

/•OO 

Nn^oih) ^ / {xc + xecxp[-{s-he)/hs])n{s)ds (A.20) 



OO 



= Xc n{s) ds + Xeexp{hi,/hs) / exp{~s/hs)n{s) ds 

Jh Jh 

= XcNe{h)+x',Ne{h), 
using the abbreviations 

Nc{h) = / n{s)ds (A.21) 

Jh 

/•OO 

Ne{h) ^ / e-''/''^n{s)ds. 

Jh 

The minimum column is given by the constant volume mixing ratio, setting 

NH,0,rmn=XcNcih). (A.22) 

Hence, 

x'e = (A^H.o(/i) - Xc Ne{h)) /N,{h). (A.23) 
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Similarly, using the abbreviations 

/•OC 

KAh) = / n{s)kn,o{'^,P{s),T{s))ds (A.24) 

Jh 
/•oo 

Jh 
/■oo 

Jh 



h 
/•oo 

K,,{h) = / n2(s)e^/''%</;(i.,P(s),r(s))ds, 

the water opacity is 

VH20W = x,{K,{h)+x,K,,{h)} (A.25) 

« x,K,{h)+x'^Ke{h). 

Note that in this notation, the absorption coefficient of water includes the contributions from CIA-processes due to 
collisions with other species, 

fcH,o(^,P(s),T(s))^fcH,o(^,P(s),T(s))+ Y. H<s)Kn,o,d^.P{S),T{s)). (A.26) 

?#H20 

The self-broadening contribution for water in the thin, upper atmosphere can be neglected (jBauer et al.l . 119871 Il989[ ). 
Combining the opacity contribution from the constant volume-mixing ratio part of water, Tc^uih) = XcKc{h), with the 
opacity of all other species into 

C#H20 

we can write the zenith opacity as separated into a dry opacity, resulting from all absorption from species other than 
water, but including the constant mixing ratio contribution from water in the upper stratosphere, and a wet opacity, 
which increases in proportion with the precipitable water vapor from the exponential part of the water vapor above the 
observer, 

PWV = -^Sl^ (A^H,0(^) - Nn,0,rmn) ■ (A.28) 

TOH2O 

Hence, we obtain equ. lA. 161 with 

^'.(M^^^^SlO 1 Xe(ft). (A.29) 

PH2O Ne[h) 

A. 6.1. Alternative parameterization of the atmospheric water vapor 

Equ. IA.19I is a rather ad hoc parameterization, which approximates the water abundance altitude profile above the 
SOFIA flight altitudes. An alternative parameterization, which is actually the one used in the present implementation 
with the kalibrate task in the kosma_software package, would be as follows: 

2^H2o(s) = I |J;^;_(,_,^)/,^)^ J;>^; (A.30) 

_ i,s-h.)/h.) , ^ f 1 - e-(^-''=)/"=), h>K 

This gives, using Ne{h) and Nc{h) as defined in equ. IA.211 above 

NuMh) ^ Xce'''^''' I I e-^n{s)ds + e~T^ I (l - e'^) n{s) ds\ (A.31) 



Xc e''=/"= [N,ih) - N,ih,) + e-''^/''^N,{K)} . 



As in the case above, he can be varied to scale the water profile to the desired total column of water above the observer. 
Owing to the different parameterization, the variation of A^H2o(^) with h^ is now a more complex dependence and solving 
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for /ig is not possible analytically. The minimum column of water above the observer is now given by the contribution 
with the constant volume mixing ratio when h^ ^ h, namely 

Nii2o\h=h = ^c n(s)ds = XcN^{h), (A.32) 

Jh 

identical to the case above. Similarly, and again neglecting the self-broadening contributions of water, the water opacity 
is now given by 

r.3,o = xee''^/'*^ [K,{h) - K,[K) + e~ '''/'' ^ K,{h)'^ . (A.33) 

For small increases above the minimum, i.e. /ig = /i + e, we have 

+x,{N,ihe)~N,ih)) 

fh+t 

n{s) ds 



j-h+c 

^^g(h+e)/h, / n{s)ds- 
Jh 


- Xc 1 ' 
Jh 


x.e^l^' (\ + -^ n{h)e- 


- Xc n{h) e 


n(h) n 
Xc—, — e ■ 





e-i-O 



For the water opacity we obtain, following eau. IA.3"0l and abbreviating k{s) = k-i{^o{v,P{s),T{s)), 

^-H.oi;:::^' = xce^''+'^/''^{Kc{h)~Kc{h+e)) (A.35) 



+Kc{Nc{hc) - Kc{h)) 

2.^e(''+')/''= / n{s)k{s)ds^xc n{s)k{s)ds 



£^0 



XcC^'^' (i + Tt) n{h)k{h)e~Xcn{h)k{h). 
n{h) 



■^c , ■k{h)e\ 

so that the opacity is proportional to pwv^ with 

b,{h) = ^5i^ A:h,o(^, P{h), T{h)). (A.36) 

For larger water columns, the opacity is not necessarily proportional to pwv: although the iVg-, resp. Kc-, contribution 
in eaus. ET341 and IA.35I becomes rapidly insignificant relative to the first terms because of the exponential decrease with 
increasing he '■ the variation of the altitude distribution due to the changing point of intersection between the constant 
and exponential part, he, results in changing line profiles in the first terms. However, as long as the variation of the line 
profile with altitude is relatively slow, i.e. ki^{s) ~ k^{h), we approximately have 

Kc{h) - Kc{hc) « k,{h) {Nc{h) - Nc{hc)) , (A.37) 

so that the proportionality between opacity and water column holds again, with equ. IA.36I holding approximately. 

A. 7. Spectrometer response and gain profile 

The spectrometer detects in each frequency channel (we omit the index specifying the spectrometer channel, resp. 
intermediate frequency, vif, in the following) a count rate (counts/sec) C proportional to the radiation input plus the 
receiver noise (expressed as an equivalent brightness temperature at the receiver input. Tree), the total being 

J^sys ^ J^rec ~^ J- A- (A.ooj 

Additionally, depending on the detailed detection mechanism of the power in the spectrometer channels, the receiver 
may have a count-offset Z , i.e., a non-zero count rate at zero (IF-) input power. The detected count rate is thus 

C = -/{Ta + Tree) + Z, (A.39) 

so that the spectrometer gain 7 (gain profile, if specified as a function of frequency) determines the response of the 
receiver for any difference measurement: 

AC = C2 - Ci = j{Ta.2 - Ta,i) = 7 ATa. (A.40) 

11 
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A. 8. Gain calibration: Hot- and cold-load measurements 

For calibration measurements on the hot- or cold-load, the mixer sees 

1 j^ hot = (CxsXg -\- LriXijlhot + 
^cold cold 

\^ -^s )^sui/s\ term) ' \^ •^i)'--^i<J Us \-^ term) 7 

where we have used the abbreviation 

GgXsJu, {Thot ) -I- GiXiJu, [Thot ) 

hot — 



(A.41) 



GsXs{l + T) 



(A.42) 



With the known physical temperatures of the loads, and the known sideband gains and sideband filter coupling coefficients, 
the spectrometer gain in each spectral resolution element can be measured from the difference in count rates between 
the hot- and cold- load measurement: 



Chot ~ C, 



cold 



Chot — c, 



cold 



Ta. 



hot 



T„, 



old 



{GsXs -\- GiXi) {{Tl^ot ~ "^coid) 



(A.43) 



As discussed above, the difference between signal- and image-sideband brightness temperatures may be neglected to 
relative high precision, resulting in Tl^t = Ju{Thot ), where an obviously consistent choice for the frequency at which to 

cold coZd 

evaluate the brightnesses would he u — ^{vs + ^i) = i^LO- 



A. 9. Sky transmission calibration 

The sky transmission is derived from measurements on the blank sky (either explicit sky-measurements, or the off-source 
positions of an astronomical observation) . The measured difference between sky- and hot- load measurement gives 



{TA,sky — Ta. 



G. 



hot) 



sky 



Ghc 



(A.44) 



(note that here and in the following, the cold-load measurement and temperatures may be used instead of the hot- load). 
This is compared with the theoretical expression for the brightness difference, which, using equations IA.6I and I A. 411 gives 



[J-A.sky ^^i/iotjljneas — ^ Jamb) X 

\^ '^a.Us)^ S"^ s -^ sky .Us ^" \^ ^a.Ui )^i'^i -^sky.u-, 
Jamh i^s-^s <J Us\-^ amb) t ^iXi Ou^\J-amb)\ 

[GsXsJusiThot) + GiXiJuXThot)] 

\^ Jamb) ^^ S'^ s ^sky^UsV a.Us ) -^hot.u, 

-\-yjlXi -LskyMiK^ ^a^Ui ) ^hot^ut 

where we have used the abbreviation for the effective hot-temperature 

r^ <Ju\-^hot) Jamb <J u\-^ amb) 

-^hot^u — Z 7 ■ 

-L Jamb 

Introducing the additional approximations 



^ sky , Us 
-^ hot.Us 



, CLi 



^ sky,Ui 
^ hot,Ui 



(A.45) 



(A.46) 



(A.47) 



r = 



GjXt 



T. 



U! 



sky Mi 



T 



sky,Vs 



(^ -^ hot.Vi -^ sky,. 
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gives 

{TA,sky - TA,hot)\,neas = ^(^ ^ /a™b) Thot.i^, a^ GgXs (1 + T) X (A.48) 



^a.jy. 



l + Tuj {tg^^Jta^^J 1-a, l + ri7 

1 + r as 1 + r 



1/q 
V-t Jamb) ^sky^Us ^s^s ^ 

ia,., (i + Tl^ (W./WJ) + ^^^ (1 + rr!) 

as 

As discussed above, w ss 1 to high precision for the relevant range of frequencies and temperatures, fi ss 1 only holds 
as long as the effective sky brightness is close to a black-body function, i.e., Tsky.v ~ Cfv{Tsky) , with a frequency- 
independent effective physical sky temperature Tsky This is not necessarily the case: a narrow- line absorption in the 
signal band resulting from water at the low temperatures of the upper troposphere and a similarly narrow feature from a 
high-altitude dry species at the higher physical temperature of the upper stratosphere may have very different effective 
sky temperatures, following their definitions according to equations lA. 171 and I A. 18l 

Assuming for simplicity fi = 1 and w = 1, and considering a double-sideband receiver with equal gain in both 
sidebands and no sideband filter, i.e., F = 1, Xs = 1 and Gg = 1, we obtain the simple formula 

{TA,sky — TA,hot)\meas ~ ^(^ ^ famb) Jus{Tsky) { ^ {'ta,v, + ^a.jyj H \ ■ (A.49) 

For a single-sideband receiver, i.e., Gi = 0, and hence F = 0,Xs = 1 and Gs ~ 1, we obtain 

{TA,sky — TA,hot)\meas ~ ^(^ ^ famb) Jv^ (Tsky) S ta,iy^ H \ ■ (A.50) 

The latter can be immediately inverted to calculate the sky transmission in the signal band; in the former case, only the 
average transmission in the signal- and image-sideband can be determined. 
The same holds for the general case of equation IA.48[ which gives 

ia,.. +F..t,., = (.Cno.-Csky)h _ i_^ (1+rfl), (A.51) 

(i^ — Jamb) J^sky,Us GgXs "s 

from which the sky transmission in the signal-sideband alone can be determined only with additional knowledge on i) 
the ratio of signal- and image-sideband transmission, and ii) the effective sky brightness at the signal- and image-band 
frequency, Tsky,v^ ;, which is also needed to calculated Vt with full precision. 

Determining the transmission across the spectrometer channel by channel is not the best strategy, because the values 
derived would be very noisy due to the relatively low signal-to-noise ratio per spectrometer channel. This can be improved 
by using broader frequency averages, but at the cost of loosing the information of variations of the atmospheric absorption 
with frequency on frequency scales below the averaging interval. 

To overcome the difficulty of the unknown atmospheric transmission ratio between image- and signal-sideband, one 
has to use additional knowledge about the atmosphere by applying an atmospheric model and fitting the observed sky-hot 
brightness to the model-predicted one. 

This is implemented in the KOSMA-observing software package kosmasoftware in the task kalihrate as follows: (i) the 
first step is to identify and discard "bad" channels, e.g., at the band edges without significant IF-response and saturated 
channels etc. (ii) the observed and calibrated sky-hot difference across the reception bandwidth of the receiver is then 
fitted with the atmospheric model prediction, following eau. lA.4^ using the pwv as a free parameter; the implementation 
at present makes the simplifying assumption of w = 1 and = 1. The best-fit value is thus determined with high 
precision, making use of the full spectral information observed, i.e., with a good signal-to-noise ratio. The user can 
control whether the fit is performed as a "common" fit to all simultaneously observed spectra (several receiver pixels at 
several frequency bands) or individually for each spectrometer, (iii) The measured on-off source spectrum can then be 
calibrated on the main-beam brightness temperature scale via 

ATjnb.ssB = T^ exp {Kpwv + c^). (A.52) 

flmb 7 GsXs 

This is made separately for each spectrometer channel, thus applying the frequency dependence of the atmospheric 
transmission as specified by the atmospheric model to the data. 
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A. 10. Determination of receiver noise and system temperature 

Equs. IA.38I and IA.39I define the system temperature used for proper noise-weighting of the data when averaging, and 
the receiver temperature that needs to be optimized when tuning the receiver, both on a Ta scale, however. Comparison 
with equ. IA.41I shows that the system temperature and receiver temperature referred to an equivalent black-body with 
physical temperature Tsys.phys, resp. Trec,phys at the receiver input, either double-sideband or single-sideband, is 

T T 

^rec,phys,DSB — ^i TT^ -Lrec,phys,SSB — T; (^A.Odj 

respectively 

Ts s Ts s 

Tsys,phys,DSB = 'Fr—~Fr— Tsys^physSSB = TTZT ' (A.54) 

From the defining equ. IA.39[ we obtain for the hot- and cold-load measurement 

C'hot - Z = 7 {Tree + TA,hot) (A.55) 

Ceold — Z = J {Tree + TA,eold), 

and can derive (Y-factor method) 

-,, (-'hot ~ Z J-ree ~i' J- A.hot ,. ,.„% 

^ = -F< : = ^? 7^F- • (A-56) 



Solving for T^ee gives 

rp _ J^ A,hot ~ i 1 A,cold 

J^ree — Y -I 



(A.57) 



^ [Ct<;(1 — XsjJv^iTterm) + (-rj:(l ~ Xi)jin{Tterm)\ 
~ ^ree\^sXs + CriXi) ~ [(^s{^ ^ Xs)jv^\Tterm) + ^i(l ~ Xi)Jiy-[Tterm)\ 

SO that 

rp _ rp, _ G,{l-x,)J^^(Tt„^) + Gi(l-x,)J^^(Tt„„,) 

J^ree,phys,DSB — J rec GsXs+Gi'Xi 

lrec.,phys.,SSB - 1 reeK^ + ^ ) gI^, 



(A.58) 



Ignoring the difference in brightness between the signal- and image-sideband on the hot- and cold load, T/g^ reduces to 
the standard formula for the receiver temperature in the y-factor method, 

rpl ^ Jv{-'-hot) ~ J Ju\^eold) / . |-„n 

^rec Y ~l ' y^-oy) 

The additional term on the right-hand-side corrects for the coupling to the image-sideband termination in the sideband 
filter and reduces to with Xs = Xi = 1. 

We obtain for the system temperature (on the off-source position on the sky) 

^ sys = J-ree ^ -I- A,sky (A. 50) 

— {Tree + TA,hot) + {TA.sky — TA,hot) 
= iKee + Uot) GsXs{l + T) + {TA,sky - TA,hot) ■ 

Insertingly, we obtain from equ. IA.48I 

Tsys = GsXs{l + r) X (A.61) 

1 1-a, 1 + rf}' 



-'- ree ~^ -'- hot ^sky.v^ V-'- Jamb) I ''a,Us + i , "n 

\ a tts I + 1 

and Tsys,DSB or Tgys^ssB following equ. IA.54I 
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